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D ■ Abstract 

QQ ■ We calculate the twist three distribution f±{x, A;-*-) contributing to Cahn effect in unpolarized semi- 






inclusive deep inelastic scattering. We use light-front Hamiltonian technique and take the state to be 
a dressed quark at one loop in perturbation theory. The 'genuine twist three' contribution comes from 



p i' the quark-gluon interaction part in the operator and is explicitly calculated. f±{x,k-^) is compared 

^ '. 
^ . with fi{x,k-^). 
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I. INTRODUCTION 

Transverse momentum dependent parton distributions (TMDs) [l| have gained a lot of in- 
terest recently. In collinear hard scattering processes, for example, in deep inelastic scattering 
(DIS) the large virtuality Q^ of the hard probe (virtual photon) introduces a longitudinal di- 
rection. A plane perpendicular to that is the transverse plane. Ordinary parton distributions 
(pdfs) measured in inclusive processes like DIS do not give any information on the transverse 
momentum distributions of quarks and gluons. TMDs can be measured in processes when suf- 
ficient transverse momentum is measured in the final state; for example in semi-inclusive deep 
inelastic scattering (SIDIS) where a hadron with transverse momentum P^ is measured or in 
Drell-Yan process where the transverse momentum of the virtual photon is measured. Factor- 
ization for some processes involving TMDs have been proven at twist two at one loop order and 
argued to hold for all orders |2|, yj. The TMDs involve an operator structure which is bilocal 
both in the light cone as well as in the transverse direction; and a path ordered exponential 
of the line integral of the gauge field (gauge link) is necessary for color gauge invariance. In 
light front gauge, the link in the light cone direction becomes unity, but contribution comes 
from the part of the link at light-cone infinity involving the transverse component of the gauge 
field. It has been found that this is process dependent in general j^ and may contribute at 



leading order in l/Q [5|. However for fragmentation process the shape of the Wilson line has 
no effect on observables |6|. In SIDIS and Drell-Yan the TMDs are simply connected by a 



reversal of sign [.?!]. In more complicated processes like hadron production in hadron-hadron 
collisions, although the standard universality of TMDs does not hold, predictivity is not lost 
|8|. However, very recently in [Oj it has been shown that such 'generalized' factorization does 
not hold for all hadropro duct ion processes. 

There are 32 quark TMDs including twist two, three and four. The higher twist or subleading 
in l/Q TMDs contain one or more 'bad' light cone component of the quark field, and the 
operator involves quark-gluon interaction term. The subleading twist TMDs are important as 
they contribute in several single spin as well as azimuthal asymmetries in the kinematical range 
of present experiments. Experimental data on several of these asymmetries are now available 



10l |. Interpretation of these subleading twist asymmetries are more challenging as they involve 



several higher twist distribution and fragmentation functions. Here model calculations of these 
functions play an important role. Twist three TMDs are related to twist two TMDs and 
'genuine twist three part' through equation of motion relations JLl|. In certain models some 
other relations between the TMDs exist based on Lorentz invariance. The Lorentz invariance 
relations do not hold in QCD due to the presence of the gauge link ^] whereas the equation 



of motion relations still hold 



ll|. 



Among the twist three TMDs there are few model calculations of f±{x, k-^). This is time- 
reversal even and plays an important role in cos (ph asymmetry in unpolarized SIDIS, the 
so-called Cahn effect 13|. The unpolarized SIDIS cross section depends on the azimuthal 
angle (ph between the lepton plane and the hadron production plane, and on the transverse 
momentum of the detected hadron. The (ph dependence of the SIDIS cross section has been 
experimentally detected by the EMC collaboration IJ]. If one neglects the explicit quark gluon 
interaction terms in the distribution and fragmentation functions then this cos cph dependence 
of the cross section is given in terms of the unpolarized distribution and fragmentation functions 
at 1/Q level. This effect has been investigated in a parton model approach by introducing a 
phenomenologically motivated intrinsic fc"*" dependence [15|. /±(x, fc"*") has been calculated in a 



simple spectator model in 16|. A bag model result of /_L(a;, fc"*") has been given in 17l |. 

In this work we calculate f±{x,k^) in light front Hamiltonian approach. Instead of using 
the Feynman diagrams, we expand the state in Fock space in terms of multi-parton light-front 
wavefunctions. The partons are on-mass shell interacting objects having non-vanishing trans- 
verse momenta and thus they can be called field theoretic partons. The advantage is that these 
wave functions are Lorentz boost invariant 181], so we can truncate the Fock space expansion 
to a few particle sector in a boost invariant way. We take the state to be a dressed quark at one 
loop in QCD. The two particle light-front wave functions (LFWFs) can be calculated analyti- 
cally for a quark at one loop using the light-front Hamiltonian. Using the constraint equation in 
light-front gauge, the bad component of the fermion field, ip^~\ is eliminated. The operator has 
a mass dependent part, a k-^ dependent part and a quark-gluon interaction part. The distri- 
bution can be expressed in terms of overlaps of LFWFs. In addition to diagonal overlaps there 
are particle number changing off-diagonal overlaps. Twist three distributions gT{x),e{x) and 
/ii(x) have been investigated in this approach before 19142 1|. In the next section, we present 



details of the calculation. We end with discussions of the result. 

II. TWIST THREE DISTRIBUTION f±{x,k^) 

The transverse momentum dependent distribution f±{x,k^)is defined as 

f/(0, y) is the path ordered exponential (link) required for color gauge invariance. For transverse 
momentum dependent distributions, the bilocality in the operator is both in the longitudinal 
as well as in the transverse direction. In the light cone gauge, A~^ = 0, the gauge link in 
the longitudinal or light-cone direction becomes unity, but contribution will come from the 
transverse gauge link at light cone infinity which can not be set to unity in this gauge. It 
has been found recently that this part of the gauge link gives important contribution even at 
twist two level, in particular in the case of time- reversal odd observables ^]. However, in the 
following, we neglect the contribution from the transverse gauge link. 

We take i = 1. We have, using the light-front projection operators A^ = |7°7^; 

In light-front gauge, A^ = 0, the 'bad' component, tjj^~'> is constrained, and the equation of 



constraint is given by 



22| 



i^^-^\y) = ^(^«^ • ^^ + 9(y^ ■ ^^ + /^HV^^+Hi/); (3) 



where the operator ^ is defined as 22 1 



1 _ 1 f°° _ _ 

Q^fix ) = 4 / dy t{x -y )f{y ). (4) 

The antisymmetric step function is given by 

eix-) = --V f—ei---- (5) 

TT J U 

V denotes the principal value. Using the equation of constraint the field i})^' can be removed. 
The operator has three parts : 

O,. = z^W^O) [(«^- d^){^h' + «'(^)(«^- d-")] ^W(y); (6) 



0, = #WT(o) {a^.A^){ 



:^)«^ + «^(^)(«'-^') 



r^Ky); 



(7) 



Or. 



m^(+)"^(0)7^ 






(8) 



For the dynamical field ^'-"'"^ we use two component formalism 22j | P, a) is a proton state of 
momentum P and helicity a. The state can be expanded in Fock space in terms of multi-parton 
LFWFs. Instead of the proton we take the state to be a dressed quark. Fock space expansion 
of such a state can be written as : 



P,a) =0i6t(P,a)|O) 



'^^J v/2(27r)3A:+ J y^2i2^.^ 



'^^^'^'^^ ^2{27TrP+6'{P-h-h) 



(f)2{P, a I h, ai; fca, X2p{h, a,)a\k2, X2) \ 0). (9) 

Here a^ and b^ are bare gluon and quark creation operators respectively and 0i and 02 are the 
multiparton wave functions. We introduce Jacobi momenta Xi, qi^ such that Yli^i ~ ^ ^^"^ 
J2i1i'^ = 0- They are defined as 

kt 



Xi 



Qi 



K: CCj 1 



Also, we introduce the wave functions, 

^1 = 01, i^2{x,,qj 



^Hkt.h^): 



(10) 



'11^ 



which are independent of the total transverse momentum P-*- of the state and are boost invari- 
ant. The two particle wave function depends on the helicities of the quark a nd g luon. Using 
the eigenvalue equation for the light-cone Hamiltonian, this can be written as 



r2.A^,q^) 



xil — x) 



1 
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{q^y + m2(l - x)2 ^{1-x) v/2(2^ 



T'-xl 



q^ 



1 — X 



X 






X 



12) 



m is t 
state 



le bare mass of the quark, cti = o"2! <^2 



23|: 



^1 



1 <^s^ . Q 

1 - ^Cf log- 



0"!. ipi actually gives the normalization of the 

2 A-l-e 



27r 



dx 



l + x^ 



X 



(13) 



to order a^. Here e is a small cutoff on x. We have taken the cutoff on the transverse momenta 
to be Q, the large scale of the process. In the above expression, we have neglected sub leading 
finite pieces, /i is a small scale such that (q^)'^ > fi^ » rn^ . 

Contribution from Om is zero. O^x has contribution from single particle sector as well as 
two particle sector of the state and Og will get contribution from an overlap of a single particle 
and a two-particle light-front wave functions. The contributions from Ok^ are : 

d_y_0^ ^iP-y-.^^.--y-^ I o,. I P) = 6{1 - x)^5\k^ - P^)\ ^, r 

ciV^^^^I i^isA^.q^) \'5\k^-q^-xP^) 
■.^6{l-x)6\k^-P^)\i„\' 



Here we have summed over the helicity of the state. The Fock space expansion of the interaction 
part of the operator can be written as 



\kjy +ik;^-y 



Of^ = 3 E J (dk^) J (dk,) J [dk,] xl (4 -^^36^3)xA. e-i' 

spins 

\-rT^bl(k^)bx,ik2)ax,{ks) + T^r^blik,)b,,ik2)aliks)\ (15) 



*' )XX2 



spins 

11.+..-. 



^-^kjy +ikj^-y^ 



^--hl{k,)h,,{k,)a,,{k,)e~'^'ty-^^^^-y' 



^2 + '^d. 



1 

"2 ^ '^3 



+ T::^^-^hl{k,)h,,{k,)al{k,)ei'ty'-H-y- . (16) 

fCo /Co 



Here we have used the notations idk) = "^^ '^. ^, ^ and \dk] = ^^-tttx- As we stated above, 

^ ^ 2{2-n-yVk+ L J 2(27r)"'fc+ ' 

the interaction part of the operator gives iplip2 and ip2'^i type terms. Contribution from Og is 
given by : 



J 4(27r)' \ \ 9\ I 27r2 ^7 ^ 



4(27r)' X I </! / 27r2 ^y P+ [(g±)^ + ^2(1 _ 3^)2] 

^^52(^±_^±_^p±)^ (17) 

a;(l — x) 



Here we have used explicit form of ip2s a(^? 1^)- 
In the frame where P^ = 0, one has 



xf±{x,k-^) 



as /{k^f{l + x^)/{l~x)+m^{l-xf 



\2 „2 I ^2/1 ^\2 



27r2 V p±)2 ^ ^2(1 _ ^)2]^ (1 - x)[{k^Y + m2(l - x)2] 



as [ik^)x^ + m\l - x)\x - 2)x] ,^^, 
27^2 {1 - x)[{k^f + m\l ~ xY] 

The twist two unpolarized distribution fi{x, k-^) can be calculated using the definition 



H^^k^) = I r£_ij^e^^"^--e-''="-^"(p | mu{o,yh^Hy~ ,y^) I P) L+=o- (19) 



The operator neglecting the gauge link is of the form 2ijj^~^'Y0)w {y~ ^ y^)- ^oi a dressed quark 
state, one gets 



h{x,k^) = 6{l-x)6\k^-P^)\ij,f 

«^^ /.2X.2.,^ ^ „K[(?Y(^)+rn^(l-^)^] 



;Cf I Sq^5\k^-q^-xP' 



27r2 y [(g±)2_^^2(l_2;)2]2 

In the frame P-*" = one gets 



(20) 



fii^,k^) = ^,C / ^^-^ ^ ^; (21) 

neglecting the single particle contribution at x = 1 and k^ = 0. The above result agrees 
with [24]. Note that in order to get the correct behaviour at x = 1 one has to include the 
single particle contribution and the normalization of the state Eq. (iT3l) . Comparing we see the 
equation of motion relation 

xf± = xh + h (22) 

is satisfied, with /^ is the genuine twist three quark-gluon interaction part which in our calcu- 
lation, comes from Og. 

In the TMDs we did not use the large /c"*" approximation. However in the limit of large fc"*". 



the twist three distribution has -4tt behaviour as shown in 25 1 
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(a) (b) 

FIG. 1: (Color online) Plots of (a) fi{x,k ) and (b) xf±(x,k ) vs x and k 
m = 0.3 GeV. k^ is in GeV. 



k 1. We have taken 



In Fig. 1 we have plotted fi{x,k^) and xf±{x,k^) as functions of x and k^. Substantial 
difference is observed in relatively lower k^ region, in fact xf± also becomes negative. This is 
due to the quark-gluon interaction contribution to the twist three distribution, and unlike the 
bag model [l7|. We took m = 0.3 GeV. We have divided both plots by ^. One has to be 
careful not to compare the numerical results of the dressed quark calculations with experimental 
data. However, the qualitative behaviour is interesting as unlike phenomenological models, the 
genuine twist three part comes from explicit calculation of the quark gluon interaction term. 

In the integrated distribution, there is an integration over fc-*-. The operator is bilocal only in 
minus direction. As a result, the gauge link is only in the light-cone direction and becomes unity 
in the light cone gauge. The operator can still be separated into three parts, 0^, Of.± and Og 
using the equation of constraint for ip^~\ Om, as before gives zero contribution. Contribution 
of Og is zero after k-^ integration, due to rotational symmetry. The entire contribution comes 
from Ofcx : 

1 +a;2 






EL 
p+ 






dx- 



1 — X 



(23) 



rhs is the twist two unpolarized distribution function /i(x, Q^)- This is expected as when 
integrated distributions are concerned, the transverse component of the bilocal current given 
by Eq. ([2]) with bilocality only in the minus direction, has the same parton interpretation as 
the plus component. Note that we get nonzero result only when P^ is nonzero 23 1. 



III. DISCUSSION 

In this paper, we calculate the twist three distribution f±(x, k-^) in light-front Hamiltonian 
approach. This distribution is known to play an important role in the observed Cahn effect in 
unpolarized SIDIS. Instead of a proton state we take the state to be a dressed quark at one loop 
in QCD. The advantage is that the higher Fock space component (two particle) LFWF can be 
calculated analytically. These play an important role in the higher twist distributions. The 
partons, that is, the quarks and gluons have non-zero transverse momenta and they interact. 
The transverse momentum dependence of the two-particle LFWF is obtained by solving the 
eigenvalue equation of the light-front Hamiltonian. At 0{as) this calculation is exact. However, 
we neglect the contribution from the gauge link at light-cone infinity. The operator has three 
parts, an intrinsic transverse momentum dependent term, a mass term and a 'genuine twist 
three' quark-gluon interaction term. Contribution from each of these terms are calculated 
using overlaps of LFWFs. The equation of motion relation connecting f±{x,k-'-) to the twist 
two unpolarized distribution fi{x,k^) and a quark-gluon interaction part is shown to hold. 
x/_l(x. A;-*") differs substantially in qualitative behaviour from fi{x,k^) in low k^ region. The 
last part vanishes when integrated ove k^ and one gets the same information as in fi{x,Q'^). 
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